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ANGLE BASED COMPUTATION WITH CUT-OFF ALGORITHM FOR MOVING TARGET 



Performance evaluations of a maneuvering target require Radars, Electro -optical tracking systems and telemetry 
systems to track the target. Angle based tracking system uses only Azimuth & Elevation which may be generated either from 
Electro-optical systems or any other antenna based tracking systems which provide Azimuth & Elevation. In such case 
minimum of any two systems are required for position calculation. For calculation of position, if only two such systems are 
considered then computation is simple. Our proposed algorithm considers two or more number of tracking systems. For such 
angle based tracking system, even if the measured azimuth and elevation are within the error limit of the instrument, still the 
problem of deviation of computed position from the true position occurs as the target moves far away from the systems. 
This is because of low apex angle generated by tracking systems. In this paper, algorithm to compute the position of the target 
is presented. Further an algorithm of Cut-off has been presented based on the computed apex angle in real time environment 
which is used efficiently to eliminate the problem of deviation. Radar measures Azimuth, Elevation & Range information and 
Telemetry system uses Inertial Navigation System for providing more accurate position information at far distance of the 
target. The performances of algorithms are analyzed based on various real time track data which are compared with Radar and 
Telemetry systems. 
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I. INTRODUCTION 

For an airborne target, position calculation is essential with respect to a particular reference point. This reference 
point may be the starting point of the target or any point towards which the target is moving. Radar based tracking system 
provides Azimuth, Elevation and Range information of the target so that a single Radar can provide position information of the 
target with respect to its own position. Then a suitable co-ordinate transformation method is used to compute the position with 
respect to the reference point. In case of angle based tracking system, a single system can provide only azimuth and elevation 
information which is insufficient to find the position of the target. In such case, minimum two systems are required for 
position calculation [1] [2]. Section II presents the algorithm for computing the position of the target for such angle based 
tracking systems. Section III deals with the analysis of positional deviation with the other sensors data such as Radar and 
Telemetry. In section IV the algorithm for apex angle computation and cut-off is presented along with the performance 
analysis of proposed apex based cut off algorithm. 
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II. COMPUTATION ALGORITHM 




Figure 1: Axis Convention ECEF and ENV 



Two frames of reference are used as in Figure 1 . 

• ECEF (Earth Centered Earth Fixed) axis with respect to centre of the earth and 

• Local i.e ENV (East(x) North(y) Vertical (z) axis with respect to any tracking systems or point of reference consider 
N no of tracking systems each of which provide two parameters 

• Azimuth (AZ): The angle measured w.r.t. north horizontally in clockwise direction (0 deg to 360 deg). 

• Elevation (EL): Elevation is the angle measured vertically upward (+ve) & downward (-ve) w.r.t. Horizontal plane 
(XY Plane). 




Figure 2: Position of a Tracking Station with its Line of Sight OL 

[n Figure2 , x^y^z^ is the position of ith tracking station w.r.t. reference point's local axis which can be obtained from 
the following steps. 

The coordinate of any point above the earth w.r.t. ECEF axis is [3] 

X ece f= (y+h)cosA,cos|ii 

Y eC ef= (y +h)cosXsinju 

Z ece f= {(1-e 2 ) y +h}sinA, 

Where X, ju & h are geodetic latitude, longitude and height respectively of the point above the earth and 
y=a/(l-e 2 sin 2 A,) 1/2 is the effective radius of the earth with the WGS-84 model [4] having the following values 

a=6378 137.0 ('a' is major axis) 
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e 2 =f (2-f) ('e' is eccentricity) 
f=l/298.257223563 (T is flattening) 

Using the above equation we can obtain the coordinate of ith tracking station: (X stllj Y stn ,Z stn ) and reference point (X re f 5 
Y ref , Z ref ) both w.r.t. ECEF axis of reference. 

Thus the coordinate of ith tracking station ( x^y^z^) wrt reference point's local axis is 







" — sin\x 


COS[X 


- 












— siviXcasy. 


— sinXsiny. 


cos X 






(i) 






. cosXcos[i 


cosXsin[x 


sinX . 









Matrix-1 

Where matrix-1 is the ECEF to Local alignment matrix and X, ju are the geodetic latitude and longitude of the 
reference point [5]. 

hi Figure 2, x m ,y m ,z m is the unknown position of target M, w.r.t. reference point's local axis which is to be found 
out from following steps . OL is the line of sight with which ith tracking station looking to target which may not exactly point 
towards the target M due to it's measurement error. 

Let I J, m[&nl be the direction cosines of line OPi t where these direction cosines are w.r.t. local axis of ith tracking 
station at point O which is computed from its measured Azimuth(AZ) and Elevation(EL) values as follows. 

I^cosEljsinAZj 
mpcosEi fCosA^ ■ 
n^=sinEf [ 

Now new direction cosines of the same line OPi parallel to local axis of common reference point is 
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Matrix-2 Matrix-3 

where Matrix-3 is Local to ECEF alignment matrix of ith tracking station in which case X and u are the Geodetic 
Latitude and Longitude of ith tracking station respectively and Matrix-2 is the ECEF to Local alignment matrix of reference 
point in which case X and u are the Geodetic Latitude and Longitude of reference point respectively [4] . 

Projection of OM on OL is Oft (Figure 2) 

OPi = (x m - Xi) li + (y m - yO mi + (z m - Zi ) n { 

NOW MPi 2 = OM 2 - OP, 2 = (x m - x0 2 + (y m - yi ) 2 + (z m - Zi ) 2 - {(x m - Xi ) l i+ (y m - Yi ) m i+ (z m - z^n,} 2 
Sum of the square of perpendicular distances contributed by total N tracking stations is: 

S = £ J ?Li(^L) 2 = lf =1 (C^ -Xi?-ty m -Yd 7 - (z m -^) 2 >-Sr =L {(x m -Xi)li+(y m -yi)mi+(z m -Zi)ni} 2 (2) 



144 



R. K. Patra, S. Kaity, J. S. Bharti & J. Ekka 



35 35 35 35 

The above sum S is minimum when - — =0, =0 and - — =0[6]. Thus from above equation (2) - — = gives the 

following equation. 

x m (N - EfU if) - ;v m SfLi h^i - *m EfU ' e«i = S?Li^(i - 1 1 ) - ^i^yt *f™f - ^ 



Similarly = gives the following equation (4) 
3 -> r m 

35 

and - — = gives the following equation 

Writing equations 3, 4 & 5 in matrix form, we get the following matrix equation. 



'N - Sf =1 If -Sf =1 iim ( 
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Thus solving above matrix equation gives the unknown position x m ^ y m ^ & z m of the target w.r.t. reference point. 

III. DEVIATION ANALYSIS 

It is observed that in case of long range tracking, due to small measurement error in both azimuth and elevation of 
each tracking station, deviation in position measurement occurs as computed by above algorithm. This deviation is very small 
and it is within the error limit if the target distance is comparatively less with respect to the measuring stations. As the distance 
of the target increases from tracking stations, this deviation increases due to measurement error [7]. In some other cases also, 
where the target is not far away from the measuring stations, deviation is also observed in some particular pattern of trajectory 
path of the target. In both cases the main reason of deviation is due to the geographical position of the flight vehicle with 
respect to measuring stations which is finally attributed to the low apex angle. The deviations as observed for both of the 
above cases will be explained subsequently. 
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Figure 3: Location of Six Different Tracking Station in XY Plane 
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GR Vs ALT:Ground Range Vs Altitude 




GR in km 

Figure 4: Example of Unacceptable GR vs ALT Plot with Deviation before Cutoff 

T Vs Elevation 
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Figure 5: Time vs Elevation Plot of Six Tracking Stations 

T Vs Azimuth 




T in Sec 

Figure 6: Time vs Azimuth Plot of Six Tracking Stations 

Figure 3 represents site locations of six tracking stations given in XY plane w.r.t. reference point which are located 
within 10 km from the reference point. The target moved from the reference point attending height of more than 40 km and 
ground range more than 145 km w.r.t. reference point whose track is given by radar, telemetry & fused track generated by six 
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tracking stations using the above computational algorithm in sec -II. The track result in Ground Range verses Altitude is shown 
in Figure 4. The elevation and azimuth values received from each of the tracking station for the corresponding track are shown 
in Figure 5 and Figure 6 respectively along the time axis in which few gaps in the plot represent track loss by the stations. In 
this case when the target is moving away from all tracking stations, it is natural that all azimuth and elevations will appear to 
converge and hence the apex angle between any two stations will gradually decrease which will cause deviation of the 
computed position as compared to radar and telemetry for small error in measurement in each tracking station [8] [9]. 
This deviation is clear from the Figure 4. 
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Figure 7: Example of Unacceptable X vs Y Plot with Deviation before Cutoff 

Figure 7 shows another example, in which the track result of target is shown in xy plane. In this case the target moved 
around the four tracking stations and reference point where the height of the target was between 400 meter and 2.5 km. 
The new positions of four tracking stations along with reference point is also shown in Figure 7 The elevation and azimuth 
values received from each of the tracking stations are shown in Figure 8 and Figure 9 respectively along the time axis. 
It is seen in Figure 7 that the fused data deviated from radar and telemetry data in particular region where the apex angle due to 
four tracking stations appears to be less. 
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Figure 8: Time vs Elevation Plot of Four Tracking Stations 
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T Vs Azimuth 
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Figure 9: Time vs Azimuth Plot of Four Tracking Stations 

IV. ALGORITHM FOR APEX ANGLE COMPUTATION AND CUT-OFF 

Figure 1 shows positions of any two tracking stations and target in space where is the apex angle between the two 
tracking stations. To avoid the deviations of the track result as shown in Figure 4 and Figure 7, it is necessary to apply cut-off 
to the computed data (Fused Data) based on apex angle between pair of tracking stations. Following steps represent the 
algorithm for apex angle computation and cut-off. 

Stepl: Prior to real time track, compute the square of distances between stations and put in the Matrix-4 
(assumed size: 6x6) where position of each station is known w.r.t. reference point as computed in equation (1) of section -II 
based on Lat, Long & Alt values of each station. Let Djjj be the square of distance from i t jy to JtJfe station, where 

1 < i < N & i + 1 < j < N, N=total no of stations. 

Since the measurement in this case is distance, Djl becomes redundant once Dq is computed and hence the Matrix-4 
becomes upper triangular as shown. 
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Figure 10: Apex Angle 6 between two Tracking Stations 
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Step 2: In real time, get the position of target X m , Y^.. from the existing algorithm mentioned in section II 

Step 3: For N tracking stations, get at most N values of distances and its square between target & tracking station if 
each station is able to track the target at the moment, otherwise make these values zero. Let X m i , X^ be these values for i'th 
station. 

Step 4: Apply the following formula for apex angle between two valid i'th & j'th tracking station using the law of 
cosines for triangle, where 9 is the apex angle. Cos B = mL " m ^ — u - 

Where l<i<N&i-\-l<j<N& & is apex angle between i t k station and station. 

Step 5: Quit from the above loop (step 2 to 4) if 8 D < apex < 172° between any two stations and consider the data 
(J^ , Y mJ Z m ) to be valid, otherwise continue the loop to search for any combinations of two stations giving apex such that 
8° < apex < 172°. If all the combinations finish without this condition then computed position of target to be considered as 
invalid due to low apex angle. 

The complexity of the above cutoff algorithm is O (n 2 ). 

The above apex based cut-off algorithm was applied to the track data of six tracking stations whose fused data was 
represented in Figure 4 and the result obtained is shown as in Figure 1 1 , in which the rest of the fused data matched with the 
telemetry and radar data. Similarly Figure 12 was obtained when the same algorithm was applied to the data of four tracking 
stations whose fused data was represented in Figure 7. 

Thus in both the cases the deviated data was removed in real time by the above apex based cutoff algorithm, resulting 
the matching of fused data with radar and telemetry in real time environment. 




Figure 11: Acceptable GR vs ALT Plot after Cutoff 
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Figure 12: Acceptable X vs Y Plot after Cutoff 
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V. CONCLUSIONS 

Algorithm proposed in section-II can compute accurate target positions throughout the complete trajectory path of the 
target for the theoretical accurate angular measurement of each tracking station. But practically this is not possible due to 
inherent error in angular measurement of each tracking station, and hence the deviation was more prominent outside the apex 
wind:w. Cut off implementation based on apex requires computation of apex angle by considering pairs of tracking stations 
which needs efficient algorithm for processing in real time from various possible combinations. The analysis was done for two 
cases such that in first case the target moved away from the reference point and tracking station with continuously increasing 
height. In this case deviation resulted in the ground range verses height plot as in Fig. 4. In the second case target moved 
around the reference point and tracking station in which case deviation resulted in XY plane as shown in Figure 7. In both the 
cases the above cutoff algorithm could able to remove the deviated data there by retaining the track data that was comparable 
up to the accuracy of 98.9% with ther sensors data such as radars and telemetry. Apex window: 8 D < apex < 172° was 
fixed based on track results of various type of trajectories followed by target where as in each case it could give satisfactory 
result. Both lower limit (8°) and upper limit (172°) ensure any of the two tracking station's line of sight do not appear to 
become parallel, in which case there will be large deviation of computed position for small error in measurement. 
However further analysis is required to achieve the above apex window based on performance of the particular tracking 
sensors and position accuracy needed in a particular application. 
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